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Abstract. We employ monomer-resolved computer simulations of model dendrimer 
molecules, to examine the significance of many-body effects in concentrated solutions 
of the same. In particular, we measure the radial distribution functions and the 
scattering functions between the centres of mass of the dissolved dendrimers at 
various concentrations, reaching values that slightly exceed the overlap density of 
the macromolecules. We analyse the role played by many-body effective interactions 
by comparing the structural data to those obtained by applying exclusively the 
previously obtained two-body effective interactions between the dendrimers [Gotze I O, 
Harreis H M and Likos C N 2004 J. Chem. Phys. 120 7761]. We find that the 
effects of the many-body forces are small in general and they become weaker as 
the dendrimer flexibility increases. Moreover, we test the validity of the oft-used 
factorisation approximation to the total scattering intensity into a product of the 
form- and the scattering factors, finding a breakdown of this factorisation at high 
concentrations. 



1. Introduction 

Soft matter systems are characterised by the simultaneous existence of two intrinsic 
structural length scales: the omnipresent atomic or microscopic scale that is associated 
with the solvent molecules and the monomers of the dissolved polymers (if any) and 
the mesoscopic scale that characterises the dissolved macromolecular aggregates as a 
whole. Depending on the physical system under consideration, the latter typically covers 
the range between several nanometers and micrometers, spanning thereby three orders 
of magnitude. The former is rather located in the domain of a few Angstrom. In 
attempting to bridge the scales all the way from the microscopic to the macroscopic 
ones, it has been proven very useful to eliminate the atomic degrees of freedom from 
view, by performing a statistical-mechanical trace over their degrees of freedom and 
constructing thereby an effective Hamiltonian that involves the mesoscopic degrees 
of freedom only pp. Although the effective Hamiltonian Tiefr greatly facilitates the 
transition to the macroscopic scales, both its construction and its interpretation have 
to be treated with care: indeed, the effective potential energy function that involves the 

X This paper is dedicated to Professor Lothar Schafer on the occasion of his 60th birthday. 
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mesoscopic degrees of freedom which appear in Tics are not true interaction potentials 
in the sense of Hamiltonian Mechanics but rather a constrained free energy which arises 
by the thermodynamic trace of the microscopic ones. 

There is a number of subtleties associated with the effective potential energy 
function that have to be taken into account when a coarse-grained statistical mechanical 
treatment of a soft matter system is employed. Two of them are particularly relevant in 
the context of calculating thermodynamic quantities and tracing out phase diagrams. 
First, the potential energy cannot, in general, be written as a sum of pair interactions: § 
the process of eliminating the microscopic degrees of freedom inadvertently generates 
higher-order, many-body potentials [21 El E] ■ Truncating the effective potential energy 
function at the pair-level constitutes the pair potential approximation, whose validity 
is not a priori guaranteed and has to be explicitly checked. And secondly, that the 
contributions to the potential energy are in general state-dependent, the most prominent 
example of the latter being the Debye-Hiickel effective pair potential that has been 
extensively employed to model charge-stabilised colloidal suspensions under certain 
physical conditions [5 . Sometimes the state-dependence of an effective pair potential 
hides precisely the effect of many-body forces and then particular care has to be taken in 
the ways in which the pair potential is employed, so as to avoid blatant thermodynamic 
inconsistencies [HI 13 El E] • 

Many-body potentials are already encountered in the realm of atomic systems, the 
Axilrod- Teller interaction ^U] being a characteristic example that has been shown to 
be relevant for the description of high-precision measurements of the structure factor 
of rare gases JT]. A formal decomposition of the effective potential energy function 
between the particles of one kind in a binary mixture in which the particles of the other 
kind are traced out has been given in Refs. |3] and |3]. Unfortunately, the treatment 
there applies only to mixtures for which the number densities of the two components 
can be varied at will, e.g., colloid-polymer or hard-sphere mixtures. It is not applicable 
to two broad categories of soft matter systems, namely charged mixtures and solutions 
of polymers of arbitrary architecture. In the former case, the number densities of the 
two components are constrained by the electroneutrality condition. In the latter, where 
one specific monomer ^21 or the centre of mass of the molecule ^1 El QEl El E| 
are chosen as effective, mesoscopic coordinates, the total number of monomers and the 
number of effective particles are coupled to each other through the constraint of keeping 
the number of monomers per macromolecule fixed. 

In charge-stabilised colloidal suspensions, three-body forces are generated by 
nonlinear counterion screening. Their effects have been examined by density functional 
theory and simulations jTHI as well as by numerical solution of the nonlinear-Poisson 
Boltzmann equation [201 1^ • It has been found that the three-body forces in this case 

§ An important exception, however, is the depletion attraction in colloid-polymer mixtures described by 
the idealised Asakura-Oosawa model. In this case, all n-th order polymer-mediated effective interactions 
between colloids vanish identically for n > 3 if the polymer-to-colloid size ratio does not exceed 
2V3/3 - 1. See Ref. for details. 
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are attractive P^OUl l^. a result confirmed by direct experimental measurements using 
optical tweezers I21j. As far as polymeric systems are concerned, the triplet forces 
in star polymer solutions have been analysed by theory and simulations in Ref. |22j . 
where it was found that they play a minor role for concentrations vastly exceeding the 
overlap density. For linear chains, on the other hand, the many-body forces appear 
to have a more pronounced effect, as witnessed by the considerable state-dependence 
of the effective pair potential that reproduces the correlation functions of concentrated 
polymer solutions ^3 ^] . The general functional form of the centre-of-mass effective 
interaction between polymer chains was found to preserve its Gaussian form, its strength 
and range being nevertheless modified within a range of ~ 10% of their original values, 
due to many-body effects [TTl HH] . 

Another polymeric system that serves as a prototype for a tunable colloidal system 
that displays a Gaussian, soft effective pair interaction is that of a solution of dendritic 
macromolecules, or dendrimers for simplicity It has been recently shown that a 

Gaussian effective pair potential can describe extremely well the scattering intensities 
obtained experimentally from concentrated dendrimer solutions |211I2H1- The Gaussian 
pair interaction has also been explicitly measured in recent computer simulations that 
employed two different coarse-grained models for the microscopic, monomer- monomer 
interactions |2S| • Nevertheless, in the approach of Ref. j2E] only two dendritic molecules 
were simulated, hence no information about many-body forces was gained. In the 
present work, we address the issue of the magnitude and importance of many-body 
effective interaction potentials in concentrated dendrimer solutions. We do not attempt 
to derive an explicit decomposition of the potential energy function into n-body terms, 
n = 2,3,4, .. .; this would require separate simulations of just n dendrimers. Instead, 
we explicitly simulate a large number of interacting dendrimers at the microscopic level 
simultaneously. We measure thereby the pair correlation functions in the concentrated 
system directly and we compare the result with the one obtained by simulating the 
same number of dendrimers as effective entities interacting exclusively by means of pair 
potentials. The discrepancies in the results from the two approaches for the correlation 
functions yield then information regarding the importance of the many-body forces of 
all orders. We find that the many-body effects are of minor importance, especially for 
flexible dendrimers. 

The rest of the paper is organised as follows. In section |21 we present our model 
and the simulation details. In section IHl we present our results for the correlation 
functions derived by the two approaches mentioned above and we discuss the magnitude 
and origin of their discrepancies. In section H] we turn our attention to the issue 
of the interpretation of the total scattering intensities from concentrated dendrimer 
solutions, and in particular to the question of the validity of the so-called factorisation 
approximation of the latter as the product of the form- and the structure factor, 
discussing the limits of applicability of such an approach. Finally, in section El we 
summarise and conclude. 
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2. The model and simulation details 

In this work, we focus exclusively on dendrimers of the fourth generation (G4). We 
model the macromolecules at the monomer-level using a simplified model that pictures 
every monomer as a hard sphere of diameter a. The bonding between the connected 
monomers is modeled by flexible threads of maximum extension ct(1 + 5). In detail, the 
potential between disconnected monomers is given by 

V„s(r) = |~ '-I'D m 

1 for r/cr > 1 

whereas bonded monomers interact via the potential 

for r/a < 1 




Kond(r) = { for l<r/(7< 1 + 5 (2) 

for r /a > 1 + 6 

The quantity 6 > serves as a control parameter of the dendrimer conformations, 
with small (5-values resulting into stiff dendrimers and large values into loose structures. 
This bead-thread model was originally introduced by Sheng et al. 27j, who kept 
a fixed value 6 = OA and examined the scaling of the radius of gyration of the 
dendrimers as a function of generation number and spacer length. The same model 
has been employed in a previous work by us, in order to systematically examine the 
evolution of the dendrimers' conformational properties with the generation number G 
[2H] • By comparing the results for various values of the parameter 5 and by performing 
a further comparison with results from a different model, we have shown that the 
conformational properties of single dendrimers are insensitive with respect to the details 
of the microscopic model. Moreover, this very simple, coarse-grained model reproduces 
the experimental scattering data for isolated dendrimers very well [2H1- As we are 
only interested in static properties, we also allow 'ghost chains', i.e., crossing of bonds 
occurring for 5 > ^/2 — l ^ 0.414 are possible. Monte Carlo simulations of this model are 
very fast, as there is no need to calculate energies; one only needs to check for overlaps, 
and additionally whether the conditions of the maximal bond extension are fulfilled. If 
one of these conditions is violated, the trial move is rejected in any case, so the search 
for further overlaps can be aborted. Furthermore, due to the very short range of the 
hard sphere interaction, neighbour lists are very effective. 

The effective pair interaction potential between the centres of mass of two G4- 
dendrimers has been determined with the help of configuration-biased Monte Carlo 
simulations of this model in Ref. [SHI- The strength of the interaction between 
dendrimers can be tuned by varying the number of generations or the parameter 6. 
Denoting by r the centre of mass separation, the 5-dependent effective pair potential 
V^g (r; 6) has been found to have a Gaussian form with small, additional corrections. In 
particular, it can be fitted by the function: 



/3 V}g\r; 5) = Co exp ( -— ) + d exp 



r 
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Table 1. The numerical values of the fit parameters of the effective pair potential 
between the centres of mass of two G4-dendrimers appearing in Eq. for two different 
values of 5. At the last column the gyration radius Rg,oo at infinite dilution is also 
shown. 



6 
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7o/o^^ 


ei 


7i/a2 


n/o 


£2 


72/0-^ 






0.1 


55.75 


9.75 


5.0 


0.9 


2.5 


0.1 


1.5 


7.2 


2.665 


2.0 


11.35 


33.0 


0.8 


10.0 


3.7 


0.0 






4.939 



where (3 = {ksT)"^ with Boltzmann's constant k-B and the absolute temperature T; 
the numerical values of the various fit parameters, depending on the choice of 6, are 
given in Table ^ Note that the precise values of the fit parameters are slightly different 
than those given in Ref . |2E] , since there we employed a more constrained fit by setting 
7o = 4i?^ ^/3, with the radius of gyration Rg^oo of the dendrimers at infinite dilution, 
and €2 = 0. The gyration radius is also shown at the last column of Table E Here, we 
considered G4-dendrimers with two different values, 6 = 0.1 and 6 = 2.0, representing 
the two extreme cases studied in Ref. f26j . 

Let p = N/Q be the number density of a sample containing N dendrimers enclosed 
in the volume Q. The definition of the overlap density of a dendrimer solution requires 
some care, as it is not a sharply defined quantity. Previous simulation studies with this 
system j2H] have revealed that the monomer density profiles around the dendrimer's 
centre of mass decay to zero at a distance rc = 1.5Rg^oo- Motivated by this fact, we 
envision every dendrimer as a 'soft sphere' of radius r^ and define the overlap density 
through the relation: || 

YP*^c = 1- (4) 

Moreover, we introduce the diameter of gyration at infinite dilution, r = 2Rg^oo, as the 
characteristic mesoscopic length scale to be used to introduce a dimensionless expression 
for the number density, pr^. In these terms, the overlap density of Eq. above is given 
by p*r^ = 0.566. The highest density in the simulation was PmaxT'^ = 0.605, slightly 
exceeding the overlap value, since pmax = l-07p*. 

For both values of 6, ten different concentrations were simulated, in particular at 
the densities p/pmax = 0.1, 0.2, . . . , 1.0. Periodic boundary conditions were employed 
throughout. At all densities, systems of 500 dendrimers were simulated, whereby each 
dendrimer consists of = 62 monomers, and the size of the simulation box was changed 
in order to modify the dendrimer number density. The minimum box length was 
Lmin = 9.384 r, yielding a system with the density pmax- The equilibration criterion 
for the system at hand requires some care, as there are is no internal energy in the 
microscopic model, since all interactions are either zero or infinity. We therefore took 

|j In the literature, there are alternative definitions. For polymer chains, for instance, the definition 
^p^Rl^l was used in Ref. [H]. 
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advantage of the fact that the effective, pair interaction V^^ (r; 5) between the centres 
of mass is known and given by Eq. (jH} with the parameters given in Tabled Hence, we 
chose to monitor the total effective pair potential energy f/*^^) (A^; 5) given by 

^ N N 

f/^^HAr;5) = i5^5^T/2)(|r.-r,|;5), (5) 

i = l jyti 

where Vij denotes the position of the i,j-th centre of mass. 

Two different starting configurations were tried. In the first one, the centres of mass 
dendrimers possessing identical microscopic conformations were placed at the vertices of 
a fcc-lattice, which was achieved without violation of the excluded volume conditions. 
This procedure is particularly useful especially at the highest density, Pmax, where a 
random distribution of the centres of mass will result with high probability into a 
forbidden state with monomer overlaps. The system was then equilibrated, monitoring 
U^'^\N;S) described above. In the second one, the dendrimers' centres of mass were 
placed in a random arrangement. Although this procedure requires a large number of 
failed attempts before an allowed configuration is found, especially at high densities, 
such configurations are possible. Once again, we monitored the total effective pair 
potential energy during the equilibration period, finding that it converges to the same 
value as the one obtained from the fcc-initial state. In this way, sufficient equilibration 
of the system was guaranteed. Finite-size effects were checked by selectively simulating 
some systems with 256 of dendrimers, in a box having a correspondigly smaller volume, 
so that the same density is achieved, and finding agreement between the two attempts. 

For 6 = 0.1, Neqaii = 10^ MC steps were used to equilibrate the system, and about 
Nrun = 2 X 10® steps to gather statistics. Statistical averages were calculated every 
^meas = 10 000 MC stcps. For 6 = 2.0, where a much larger random displacement 
for the monomers can be used, the equilibration phase consisted of Neq^n = 10^ steps 
and statistical averages were calculated every A^'^eas = 1000 steps over a period of 
^mn = 2 X 10* steps. The quantities measured were monomer profiles around the 
centres of mass, the radial distributions functions of the latter, radii of gyration, form 
factors, structure factors from the centres of mass and total scattering intensities; all 
these quantities will be precisely defined in the sections that follow. 

In Figs.Hand|2l we show simulation snapshots of the monomer-resolved simulations 
for the lowest and the highest density for the thread length 6 = 0.1. (For clarity, in Fig. 
^we show only a section of the simulation box of the same size as in Fig. El) Although 
at Fig. ^individual dendrimer molecules can still be resolved, since the density is much 
smaller than p^,, in Fig. |21this is not any more possible. Here, p = l.OTp* and the whole 
system appears as a dense solution of monomers, in which the individual character of 
each macromolecule is lost. We will return to the implications of this fact in section |3] 

In addition, a different kind of Monte Carlo simulations was also carried out, in 
which the monomers were not explicitly resolved. Instead, the dendrimers were replaced 
entirely by their centres of mass, which were then treated as effective, soft particles 
interacting exclusively by means of the pair potential of Eq. 0. Accordingly, we call 
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Figure 1. A snapshot from the monomer resolved-simulation of dendrimers. The 
monomers are rendered as spheres of diameter a. Here, dendrimers with threads 
characterised through i5 = 0.1 at a density pr^ = 0.0605 are shown. Note that only a 
part of the simulation box is shown, which has the same size as the full box depicted 
in Fig.ia 




Figure 2. Same as Fig. ^ but at density pr^ — 0.605. Here the complete simulation 
box is shown. 

this approach an effective simulation. As all monomers have dropped out of sight in the 
effective approach, it is only possible to measure quantities pertaining to the centres of 
mass, i.e., their radial distribution functions and structure factors. Comparison of the 
results regarding these quantities that are obtained through the two different types of 
simulations yields important information by way of testing whether the pair-potential 
approximation is meaningful. 
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3. Comparison between the monomer-resolved and the effective simulations 

Each dendrimer of the fourth generation consists of = 62 monomers. Let a, /3 be 
monomer indices within a given dendrimer whereas j are integers describing the 
dendrimer molecules as whole entities. In particular, let stand for the position of 
the centre of mass of the i-th dendrimer, RJ^ denote the position vector of the a-th 
monomer in the z-th dendrimer, and stand for the same quantity but now measured 
in a coordinate system centred at r,. Obviously, it holds 

= + K. (6) 

In the monomer-resolved simulation, the following quantities were measured: The 
radial distribution function g{r) between the centres of mass, defined as 

.INN \ 
\j=l jVi / 

where (■ ■ ■) denotes a statistical average and = rj — r^. Related to this quantity is 
the structure factor S{q) that describes the correlations between the centres of mass in 
reciprocal space and it is given by 

.INN \ 

\j=i j=i I 
Note that S{q) and g{r) are related by a Fourier transformation [2^1 

S{q) = 1 + p j dV exp [-iq • r] [g{r) - 1] . (9) 

Moreover, we took advantage of the microscopic nature of the simulation to measure the 
dendrimers' form factor F{q) at every simulated density p. This quantity is expressed 
by the relation: 



i=l \a=l 13=1 I 



(10) 



Another quantity of interest is the monomer distribution around the centre of mass, 
^(m), which can again be measured at any desired overall density and is given by the 
expression: 

i=l \a=l I 

The overall size of the dendrimer is characterised by its radius of gyration Rg, which 
was measured in the simulation by calculating the quantity: 



1 ^ 



tl\ 



(12) 



In Eqs. ()10|) - ()12|) above, the summand in the sum over i is the corresponding quantity 
(form factor, density profile, and radius of gyration, respectively) of the i-th dendrimer. 
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The additional summation over i and the division by the total number of dendrimers 
corresponds to an additional average over all dendrimers. Since all macromolecules are 
equivalent, the expectation values are identical for every summand. Finally, we also 
measured the scattering function /(g) of the concentrated solution, which corresponds 
to the coherent contribution of the total scattering intensity in a SANS experiment, 
under the assumption that all monomers possess the same scattering length density 
[SO. .3T| EHl- This is given by the equation: 

^ I N N u u \ 

^(^) = ]^ ( E E E E [-iq ■ - ^'^)] ) ' (13) 

\ i=l j=l a=l f3=l I 

i.e., it is the total coherent scattering intensity from all monomers of the system. 

In the effective picture, all information regarding the monomers' degrees of freedom 
is lost, hence in the effective simulation we can only measure the corresponding radial 
distribution function gesij") and the structure factor 5*65 (g) of the centres of mass. These 
are given by Eqs. ((Tj) and (jH)) above but with the averages now performed with the 
effective Hamiltonian, i.e., 

^eff(r) = 4(>:>:<5(r-r,,)) , (14) 




Wo 



and 




5'efr(g) = ^( ^^exp[-iq-(r, -r^)]^ . (15) 

The effective Hamiltonian Tiefr involves the momenta pj and positions rj of the centres 
of mass only and contains exclusively pair interactions, i.e., 

N 2 ^ N N 

^ef^ = E 1^ + 2 E E ^eS?^(ir^ - r.i; (16) 

1=1 i=l j^i 

where m is the dendrimers' mass, which is irrelevant as far as static quantities of the 
system are concerned. A particular property of the effective description of a complex 
system is that it leaves all correlation functions between the coarse-grained degrees 
of freedom invariant provided that the mapping into the effective system is exact pp. 
In other words, if the effective Hamiltonian contains the contributions to the effective 
potential at all orders, it makes no difference whether one calculates quantities such as 
g{r) or S{q) in the original, microscopic description or in the coarse-grained one. As our 
effective Hamiltonian Hefr is truncated at the pair level, the deviations between g{r) and 
Qesij) or, equivalently, between S{q) and S'eff(g) will be a measure of the importance of 
the neglected many-body terms in Eq. ()16|1 . 

Representative results comparing between the two approaches are shown in Fig. El 
pertaining to the dendrimers with 5 = 0.1 and in Fig.lH which refers to dendrimers with 
5 = 2.0. The length scale used in this plot is the zero-density gyration radius of the 
dendrimers, -Rg,oo- For clarity, only the results only for three different densities obtained 
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Figure 3. Comparison between the results from the monomer-resolved and the 
effective simulation of concentrated dendrimers with maximal thread length (5 = 0.1 of 
the bonds. The three different densities are p = O.lpmax, O.Spmax and pmax, as indicated 
on the plots, with PmaxT^ — 0.605. Results are shown for (a) the radial distribution 
function g{r) and (b) the structure factor S{q) of the centre of mass coordinates. 



from the monomer resolved simulations are compared to those from the effective ones. At 
sufficiently low densities, p = O.lpmax, the results from the two types of simulations are 
indistinguishable, hence the pair-potential approximation is an excellent one and many- 
body forces seem to play no role there; they can be thus safely ignored. Deviations 
between the two descriptions arise nevertheless as the overall concentration of the 
solution grows. Referring to Fig.lHfa), we see that for the 5 = 0.1-dendrimers, which have 
a rather high internal monomer density, the deviations are already visible (but small) at a 
density p = O.Spmax and they become more pronounced at the highest simulated density, 
p = Pmax- The true radial distribution function g{r) between the centres of mass shows 
a more pronounced coordination than the effective one, gc,g{r), and this effect is also 
reflected in the corresponding structure factors. The peak height of S{q) is higher than 
the one of Ses{q), pointing to the fact that the zero-density pair potential underestimates 
somehow the strength of the repulsions between the dendrimers' centres of mass. The 
relative deviation between the two descriptions as far as the peak height is concerned are 
at the highest density about 6%. Much more drastic is the discrepancy of the S{q — > 0) 
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Figure 4. Same as Fig. |3but for thread length 5 = 2.0. 



limit, for which S{q — 0) = 0.018 whereas ScE{q 0) = 0.033. Given the fact that the 
S{q = 0)-value is proportional to the osmotic isothermal compressibility of the solution, 
employing the effective picture can lead here to serious errors in the calculation of the 
thermodynamics of the system. Two integrations of the inverse compressibility are 
needed in order to obtain the Helmholtz free energy of the solution, hence errors at all 
lower densities accumulate in performing such an integration and they can lead to a 
serious underestimation of the free energy if the effective picture is employed. 

The agreement between the microscopic and the coarse-grained approaches is a lot 
better for the case of the 6 = 2.0-dendrimers, which possess a much lower internal 
monomer density than their 6 = 0.1-counterparts. Indeed, as can be seen in Fig. 
I3ja), the radial distribution functions g{r) and Qcsif') barely show any difference, all 
the way up to the maximum density pmax- Similar to the case 6 = 0.1, g{r) shows a 
slightly more pronounced coordination than Qesi^), the difference between the two is 
nevertheless extremely small. The same holds for the structure factors S{q) and Scs{q), 
shown in Fig.^fb). Here, even the discrepancy in the compressibility is very small, with 
S{q — > 0) = 0.132 and Ses{q ^ 0) = 0.138 at p = pmax- For dendrimers with a higher 
degree of internal freedom, the pair potential approximation holds all the way up to 
the overlap concentration. In this respect, it is very satisfactory that it is precisely the 
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u/g 

Figure 5. The radial monomer density profiles ^(m) [Eq. of the dendrimers 

around their centres of mass at infinite dilution [p = 0) and at the highest density 
P = Pmax = l-07p,, as indicated on the plots, (a) For model dendrimers with thread 
length (5 = 0.1 and (b) for 5 = 2.0. Note the shrinkage and growth of the profiles. 

model with the value 5 = 2.0 that has been found to accurately describe scattering data 
from real dendrimers ^Hl. 

Let us now try to to obtain some physical insight into the mechanisms that cause 
the true correlation functions to show higher ordering than the effective ones. Suppose 
that the reason lied in the increasing significance of three-body effective forces. Three- 
body potentials arise through three- dendrimer overlaps: the region of space in which 
three spherical objects simultaneously overlap is overcounted when one adds over the 
three pair interactions and it has to be subtracted anew. Given the fact that any 
overlap between repulsive monomers gives rise to a correspondingly repulsive interaction, 
together with the fact that the contribution from the triple-overlap region has to be 
subtracted, leads to the conclusion that triple forces should be attractive, as for the case 
of star polymers [221, ^^^^ self-avoiding polymer chains 02}, for which three-body 
forces have been measured explicitly j31]. Yet, an attractive contribution to the potential 
energy leads to a reduced effective pair repulsion. This is on the one hand intuitively 
clear and, on the other, it can be put in formal terms by making a density expansion 
of the density- dependent pair interaction up to linear order in density, see Eq. (10) of 
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— 8=0.1 




Figure 6. The dependence of the dendrimers' radius of gyration on the solution 
density for the two types of model macromolecules, as indicated in the legend. 

Ref. \n\- Thus, we would then obtain a weakening of the correlations and an increase 
of the osmotic compressibility, whereas in Figs. 01 and El exactly the opposite is true. In 
order to obtain the true g{r) at p = p^ax for the 6 = 0.1-dendrimers, a renormalised 
effective pair potential V^g (r; 6, p) can be employed that is more strongly repulsive than 
the original one, Vg (r; 5); as a matter of fact, we were able to reproduce g{r) at pmax 
by using (r; 5 = 0.1, pmax) — 1-2 V^efr ('"i = 0.1). A similar effect has been observed 
for polymer chains jJJ., for which the density-dependent, renormalised pair potential 
necessary to reproduce g{r) at high concentrations was found to be more repulsive than 
the one that holds at p = 0, whereas, at the same time, the correction arising from 
triplet forces alone goes in the opposite direction of weakening the pair repulsions. 

The above considerations point to the fact that the deviations between g{r) and 
gcsi'f') are a genuinely many-body effect that arises from the high concentration of the 
solution per se and cannot be attributed to three-body forces alone. In particular, 
the presence of many dendrimers surrounding a given one in the concentrated solution, 
gives rise to a deformation of the dendrimer itself. To corroborate this statement, we 
have measured the concentration-dependent monomer density profiles ^{u) around the 
dendrimers' centre of mass, given by Eq. (Illj). Results are shown in Fig. Efa) for the 
case 5 = 0.1 and in Fig. Efb) for the case 6 = 2.0. It can be seen that as a result of the 
crowding of the dendrimers at the highest concentration, the monomer profiles become 
slightly shorter in range and they grow in height; in other words, the dendrimers shrink as 
a result of the increased overall concentration, as can be also witnessed by the reduction 
of their radius of gyration shown in Fig. The molecules that effectively interact 
are stijfer at higher densities than at lower ones; their internal monomer concentration 
grows with p and as a result of this deformation, the interaction between two dendrimers 
becomes more repulsive than at zero density. 

The above claim is supported by the fact that the effect of the concentration 
on the pair interaction is much more pronounced for the dendrimers with the short 
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thread length than for those with the longer one. Although the monomer profiles for 
both dendrimer kinds grow with p, the internal monomer concentration for the stiffer 
dendrimers is much higher than the one for the softer ones. A concentration-induced 
increase of ^{u) has a much stronger effect for the effective interaction of the stiff 
dendrimers than for the soft ones, since it occurs at a scale of a^^{u) ~ 0.4 for the 
former but at a scale of a^S,{u) ~ 0.1 for the latter, see Fig. The monomer beads are 
modeled here as hard spheres. The change in the free energy of a hard-sphere fluid upon 
an increase of the local density is highly nonlinear and grows rapidly with increasing 
packing fraction, hence the effect is much more pronounced for the case 5 = 0.1 than 
for the case 5 = 2.0. 

Another way of expressing the vast discrepancy in the monomer crowding of the 
two systems is to look at the monomer packing fraction ?7m. As there are u monomers 
per dendrimer, this quantity is given by the expression 



For both types of dendrimers, u = 62 and PmaxT^ = 0.605. Yet, the ratio a/r has 
the value 0.188 for 6 = 0.1 and 0.101 for 6 = 2.0, see the last column of Table ^ 
Accordingly, at p = pmax we obtain rj^ = 0.13 for 6 = 0.1 but rj^ = 0.02 for 6 = 2.0. 
The soft dendrimers have a much lower monomer packing fraction at p* than the stiffer 
ones, a result that can be traced to the fact that their radius of gyration is larger.^ Thus, 
we conclude that the density- dependence of the pair interaction can be traced back to 
the shrinking of the dendrimers, a phenomenon that leads to increased crowding of the 
monomers in their interior. 

4. Total scattering intensities and the factorisation approximation 

In this section we turn our attention to a different question, which is however related 
to the issues discussed above, namely to the interpretation of scattering data from 
concentrated dendrimer solutions. As a first step, we consider the form factor F{q), 
defined by Eq. ()10|) . Clearly, F{q) expresses the mtramo/ecu/ar correlations between the 
monomers belonging to a certain dendrimer. When scattering from an infinitely dilute 
solution, F{q) offers the only contribution to the coherent scattering density. Since all 
the information about the monomer correlations is encoded in F{q), great experimental 
effort is devoted to the determination of this quantity. At low values of q, qRg^oo, 
the form factor delivers information about the overall size of the molecule whereas at 
higher values of the scattering wavevector, q ~ 1/a, where a is the monomer length, 
information about the monomer correlations and the fractal dimension of the object is 
hidden [H ESI ESI- 

Although F{q) is experimentally measured at the limit p — 0, the same quantity 
can be defined at any density. At arbitrary concentrations, F{q) will in general 

% This is characteristic for non-compact objects: for polymer chains, e.g., one obtains r/m ~ R;'/' at 
the overlap concentration 
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Figure 7. The form factors measured in the monomer-resolved simulations for one 
isolated dendrimer molecule [p — 0, solid line) and at the highest density [p — Pmax) 
dotted line). The model dendrimers have maximum thread length (a) (5 = 0.1 and (b) 
5 = 2.0. 



change with respect to its form at infinite dilution, due to possible deformations of 
the macromolecules. In Fig. [3 we show the form factors for the two model dendrimers 
at the lowest and at the highest simulate densities. It can be seen there that there is 
only a small change in both cases, which takes the form of a slight extension of F{q) to 
higher g-values as the concentration increases. This is consistent with the shrinkage of 
the dendrimers and the corresponding decrease of the gyration radius. Indeed, in the 
Guinier regime, qRg < 1, the form factor has a parabolic profile, F{q) = A^[l — (gi?g)^/3], 
and a reduction of Rg manifests itself as a swelling in g-space and vice versa |36j . 

Let us now turn our attention to the total coherent scattering intensity from all 
monomers, I{q), given by Eq. (fT^ . It is clear from its definition that I{q) can also 
be measured in the monomer-resolved simulation and this has been done for both 
dendrimer species, characterised by the maximum thread extensions S = 0.1 and S = 2.0. 
In attempting to model complex polymeric entities as soft colloids, it is a common 
procedure to separate the intramolecular from the intermolecular correlations and to 
write down approximations for the quantity I{q) in which the two types of correlations 
appear in a factorised fashion. Here we are going to put this approach into a test and 




Figure 8. The total coherent scattering intensity I{q) [Eq. H13|) ] from concentrated 
6 = 0.1-dendrimer solutions, compared with the result from the factorisation 
approximation, Eq. I|22|) . at different overall concentrations p. (a) p = O.lpmax and 
(b) p = O.Spmax- Results using both the form factor F{q) at the given density and its 
counterpart at infinite dilution, Fo{q) are shown for the factorisation approximation. 



figure out the limits of its validity as far as dendritic molecules are concerned. A similar 
test has been carried out by Krakoviak et al. ^0J\ who compared results from the PRISM 
model for polymers with simulations and with the factorisation ansatz. 

As a first approximate step, one assumes that the intramolecular conformations 
and centre-of-mass correlations decouple from each other. Correspondingly, Eq. (jlHj) 
takes the approximate form: 

N N V u 

^(^) = ]^ E E E E <^^P [-iq • (r. - r,)]) (exp [-iq ■ « - u^,)] > . (18) 

i=l j=l a=l f3=l 

The approximation inherent in Eq. (fTHj) above is a reasonable one for dendrimers. 
Indeed, as it has been shown in Ref. the monomer degrees of freedom are correlated 
at length scales ~ a, whereas for the overall densities p considered here, the centers 
of mass are correlated at lengths at least ~ Rg and the two are well-separated from 
each other. Hence, at the wavevector-scale gcM ~ ^/Rg which the centre-of-mass 




Figure 9. (a) Same as Figs. |Hta) and (b) but for p — p,„ax- (b) The true 
structure factor S{q) between the centres of mass at p = pmax, as obtained from the 
monomer- resolved simulations, compared with the apparent structure factor S'app((z) = 



S{q) shows structure, the dendrimers still appear as compact objects and the internal 
fluctuations can be decoupled from the intermolecular ones. The second approximation 
is now the following. Suppose that we are at sufficiently low densities, so that close 
approaches between the centres of mass of the dendrimers are very rare and they carry 
therefore a negligible statistical weight. Then, since monomers belonging to different 
dendrimers stay far apart, it is reasonable to assume that the deviations from their 
respective centres of mass are uncorrelated. In this case, one can approximately write: 

a=l /3=1 a=l (3=1 

= (19) 

where i^q is the Fourier transform of the monomer density operator ,^(u) around the 
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centre of mass of an arbitrary dendrimer:"'" 

V 

e(u) = 5^5(u-<). (20) 

Clearly, the right hand side of Eq. |T!?|l has no dependence on the dendrimer index. 
At the same time, it has been shown in Ref. [38j that the product z/~^(^q)(^_q) is an 
excellent approximation for the form factor -F(g) of the dendrimers, deviations from 
the exact expression in Eq. (fTUj) . F{cj) = {iqi-q) , appearing only at high g- values 
that are unreachable in a typical SANS experiment. The approximation inherent in Eq. 
fll9|) has been derived for monomers belonging to different dendrimers {i 7^ j) and now, 
in view of the results of Ref. j^, it can be also applied to the case i = j. Putting 
everything together, we obtain 



U V 



^ E E (^^p [-iq ■ K - > = (21) 

Eqs. (I18|) and (j21|) now yield the oft-employed factorisation approximation: 

I{q) - S{q)F{q), (22) 

whose validity will be tested in what follows. 

The assumptions that went into the derivation of Eq. (|22|) above become exact when 
the particles from which one scatters are rigid colloids [39 , in which case individual 
scattering centres are devoid of a fluctuating nature. In this context, it is important 
to note that there is an analog of the factorisation approximation that is applied 
in the theory of concentrated polymer solutions and carries the name "rigid particle 
assumption" |H71 1^ . Here, one starts from Eq. (|18p and assumes that monomer- 
monomer correlations between monomers belonging to different polymers are identical 
to the intramolecular correlations in any chain jj37|j. Under this assumption, the second 
factor on the right-hand-side of Eq. (fTHj) above takes the form: 

^ E E (-P [-iq ■ i< - ^,)] > = ^ E E (^-p [-iq ■ i< - ] > = ^(^)' (23) 

a=l (3=1 a=l 13=1 

and, in conjunction with Eq. ()18j) . the factorisation approximation of Eq. ()22j) follows 
once again. Krakoviak et al. tested the validity of Eq. (j^^ for polymer solutions, finding 
that it breaks down for high polymer densities. 

We have put the validity of Eq. into a strong test by comparing the directly 
measured total coherent scattering intensity /(g) with the product F{q)S{q), where for 
the latter quantity both factors are the ones measured in the same simulation. Results 
are shown in Figs. |H1 and IHfa) for the 6 = 0.1-dendrimers as well as in Figs. ITUl and 
ITTT a) for the 6 = 2.0- dendrimers. It can be seen that the factorisation approximation is 
valid at the lowest density shown (p = O.lpmax) but that its quality becomes poorer as 
the concentration of the solution increases. A dramatic breakdown can be seen in Fig. 

+ The quantity ^(u) defined in Eq. is simply the expectation value of the operator i^(u). 
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Figure 10. Same as Fig.|Slbut for 5 ~ 2.0-dendrimers. 

ini^a) for the more compact dendrimers, whereas the breakdown is also clear (but less 
spectacular) for the more open dendrimers, Fig. ITlT a). 

We can now trace back to the physical origins of the breakdown of the factorisation 
approximation, Eq. (j22j) . There is first of all a weak breakdown of the first assumption, 
Eq. (fTHj) . in which the centre-of-mass coordinates were decoupled from the fluctuating 
monomers. Indeed, were this approximation to be true, then the form factor F{q) would 
remain unchanged at all concentrations. This is however not the the results 

in Fig. 13 demonstrate: the dendrimers shrink as p grows. Yet, the difference between 
the infinite-dilution form factor, -Fo(g) and its counterpart at finite density, F{q), is not 
sufficient to account for the failure of the factorisation approximation. As can be seen in 
Figs. IHTb) . IWa) . fTI¥ b) and irUT a). the product S{q)F{q) is in even worse agreement with 
I{q) than the product S'(g)Fo(g). The reason for the breakdown of Eq. (j^^ lies in the 
assumption inherent in deriving the approximation of Eq. (j2^, namely that fluctuations 
between monomers belonging to different dendrimers are uncorrelated. At sufficiently 
low densities p, this is a reasonable assumption. However, in approaching the overlap 
density p=K, it does not hold any more. As monomers from different dendrimers begin to 
crowd with one another, their coordinates with respect to their centers of mass become 
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Figure 11. Same as Fig.|51but for 5 — 2.0-dcndrimcrs. 



more and more strongly correlated and Eq. (j2T|l loses its validity. In this respect, it 
is not surprising that the breakdown of Eq. (j22j) is more dramatic for the 5 = 0.1- 
dendrimers than for the 6 = 2.0-ones. In the former case, the monomer packing fraction 
is higher and the corresponding correlations between monomers belonging to different 
molecules stronger than in the latter. To put it in more pictorial terms: at the overlap 
concentration it is not any more possible to tell to which dendrimer a monomer belongs, 
see Fig. |21 A clear separation between intra- and inter-dendrimer fluctuations is not any 
more possible. 

We finally discuss the consequences of the above findings for the interpretation 
of scattering data obtained from concentrated dendrimer solutions. The validity of 
Eq. ()22j) is often taken for granted: the form factor F{q) is measured in a SANS- 
or SAXS experiment at low concentrations and extrapolated to infinite dilution to 
obtain the quantity Fo{q). Thereafter, the measured coherent scattering intensity at 
any concentration, /(g) is divided through -Fo(g), the result being interpreted as the 
structure factor of the system. In order to differentiate it from S{q), we emphasise here 
that this is only an apparent structure factor S'app(g), given by 

m 



Foiq)' 



(24) 
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In Figs.infb) and llir b) we compare the apparent structure factors for the two dendrimer 
species at the highest simulated density with the true ones. It can be seen that the 
process of applying Eq. has the effect of producing apparent structure factors that 
are everywhere lower than the true ones and they even fail to reach the asymptotic value 
unity at the range considered. 

Such structure factors from concentrated dendrimer solutions have been published 
in Refs. jH] and [l^); in which they have been correctly termed 'apparent'. It is 
important here to point out that apparent structure factors can lead to false conclusions 
regarding the validity of the pair potential approximation in mesoscopic theories of 
dendrimer solutions. Indeed, as we have explicitly shown in this work, many-body 
effective potentials play only a minor role in concentrated dendrimer solutions, therefore, 
one can obtain accurate structure factors from theory by working with a density- 
independent pair potential. If, however, these structure factors were to be compared 
with the apparent experimental quantities 5'app(g), discrepancies of the kind shown 
in Figs, in^b) and ITTT b) would show up. It would be then possible to argue that 
these discrepancies are due to the breakdown of the pair potential approximation but, 
as we have shown here, this conclusion would be unwarranted. The reason for the 
disagreement between theory and 'experiment' would, in this case, lie in the employment 
of an erroneous approximation, Eq. ()22|) . in deriving apparent structure factors from 
the experimental data. It is worth noting that Krakoviak et al. fTj] reached similar 
conclusions for the case of polymer solutions, although they did not formally introduce 
an apparent structure factor into their considerations. 

5. Summary and concluding remarks 

We have carried out extensive, monomer-resolved and effective simulations of model 
dendrimers in order to calculate correlation functions between the centres of mass of 
the macromolecules and the individual monomers themselves. By comparing the real- 
space correlation functions obtained by the two simulation approaches, we found that 
many-body effective potentials play a minor role up to the overlap density and they 
can be altogether ignored for open dendrimers with long bond lengths. Our finding for 
the scattering intensity, on the other hand, is that the factorisation approximation of 
this quantity into a form- and a structure factor loses its validity as one approaches the 
overlap concentration. Structure factors that are obtained from experimental data by 
dividing the scattering intensity through the form factor can be seriously in error. 

It appears, therefore, that the extraction of an accurate structure factor from 
concentrated dendrimer solutions is extremely difficult as one approaches the overlap 
concentration. We anticipate that this result is also valid for other 'polymeric colloids' 
such as star-shaped polymers and brushes. One strategy to circumvent this inherent 
difficulty is to use the labeling technique, in which a small, inner part of the molecule 
is protonated and the rest is deuterated in such a way that the contrast between the 
outermost part of the molecule and the solvent vanishes. In this way, only the innermost 
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part of the molecule will have contrast with the solvent and scatter coherently. Thus, 
one can reach concentrations for the whole system that exceed p*, whereas the labeled 
parts are still nonoverlapping. Such a technique was successfully applied, e.g., to star 
polymers [T^ . 
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